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It is well-known under the name of 'periodic homogenization' that, under a centering 
condition of the drift, a periodic diffusion process on R d converges, under diffusive 
rescaling, to a d-dimensional Brownian motion. Existing proofs of this result all rely 
^ • on uniform ellipticity or hypoellipticity assumptions on the diffusion. In this paper, we 

' considerably weaken these assumptions in order to allow for the diffusion coefficient 

, to even vanish on an open set. 

• As a consequence, it is no longer the case that the effective diffusivity matrix is 

CN| ' necessarily non-degenerate. It turns out that, provided that some very weak regularity 

■ conditions are met, the range of the effective diffusivity matrix can be read off the 
shape of the support of the invariant measure for the periodic diffusion. In particular, 
this gives some easily verifiable conditions for the effective diffusivity matrix to be of 

■ full rank. We also discuss the application of our results to the homogenization of a 
class of elliptic and parabolic PDEs. 
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1 Introduction 

Our goal is to study, by a probabilistic method, the limit as e — > of the solution 
u e (t, x) of an elliptic PDE in the regular bounded domain D C R d 



£ E u £ (x) + f[x,-)u £ (x) = 0, x E D, 

u E (x) = g(x), x € 3D, 

where / is bounded from above, and g is continuous, as well as the limit of u £ (t, x), 
the solution of a parabolic PDE of the form 

C s u%t, x) + f-e(-) + f(x, -)V(t,s) . 



dt \e e' e J (1.2) 

u e (0, x) = g(x), x E R d . 

In both cases, the linear operator C £ is assumed to be a second order differential oper- 
ator with rapidly oscillating coefficients given by 

1 v-^ ,x. d 2 v^fl, ,%s ,x s ~\ d 



1 , x s v r 1 1 , x s , x x 



ZJ — 1 l — l 



dxi 
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The novelty of our result lies in the fact that we allow the matrix a to degenerate (and 
even possibly to vanish) in some open subset D of R d . There is by now quite a vast 
literature concerning the homogenization of second order elliptic and parabolic PDEs 
with a possibly degenerating matrix of second order coefficients a, see among others 
HDASC92I . IBB98I . 1IBMT96I . CPPW06I . l|Par04l . But, as far as we know, in all of 
these works, either the coefficient a is allowed to degenerate in certain directions only, 
or else it may vanish on sets of Lebesgue measure zero only. It seems that our paper 
presents the first results where the matrix a is allowed to vanish on an open set. The 
main technical difficulty that we have to overcome is the lack of regularisation since 
we do not assume C e to be hypoelliptic (not even on a set of full measure). However, 
it turns out that it is possible to show nevertheless that under very weak assumptions, 
its resolvent maps C 1 into C 1 (see Lemma l2~6i l. which provides a C 1 solution to certain 
Poisson equations, and is sufficient to make an approximation argument work (see 
Lemma l3~2l >. 

Because of the high degree of degeneracy allowed by our approach, it is no longer 
obvious that the effective diffusivity A of the homogenized operator 

1 d ffi d f) 

2 ^ dxidx, f-f dxi 
i^j—i i—i 

is non-degenerate. We shall therefore also seek to characterize the image of the ho- 
mogenized diffusion matrix. It turns out that this can be done in terms of the support of 
the invariant measure of the diffusion process on the torus T d with drift b and diffusion 
matrix a. 

The paper is organized as follows. Section 2 contains our assumptions and several 
important preliminary results. Section 3 presents the homogenization result, in proba- 
bilistic terms. Section 4 contains our characterization of the image of the homogenized 
diffusion matrix A, and sections 5 and 6 present the application to elliptic and parabolic 
PDEs. Finally, Section 7 contains a few concrete example that illustrate the scope of 
the results in this paper and highlight the differences with the existing literature. 



2 Assumptions and preliminary results 

Given e > 0, x € R d , let {X* ,e } denote the solution of the SDE 

x " -*+f, H¥) + <¥)]* + gjf «(¥■) mi - <zi) 

where b, c, and a are periodic, of period one in each direction, and the process {Wt = 
(Wf, . . . , WF 1 ), t > 0} is a standard m-dimensional Brownian motion. 

Define Xf' £ = -X x f.. Then there exists a standard m-dimensional Brownian 
motion {Wt}, depending on e (but we forget that dependence since it has no incidence 
on the law of the process), such that 

m „t 

ds + J2 ^Xf s )dWl (2.2) 
j=i Jo 

In the sequel, we shall consider the solution of ( 12.21 ). as taking values in the torus T d . 
We will also consider the same equation starting from x, but without the term ec in the 



X' 



x 

F. 



b(Xf e ) + ec(X*' E ) 
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drift, namely 



X* = x- 



/ b(X*)ds + V / a 3 {X x s )dWl. 

JO 1 = 1 "'0 



(2.3) 



We denote by Jf the Jacobian of the stochastic flow associated to Xf, that is 
{Jfi t — 0}' me ^ x d-matrix valued stochastic process solving 



dJ? = Db(X?)J? dt + J2 D°j(XDJ? dWi , J$ = 1 



(2.4) 



To the SDE satisfied by the process {X®}, we associate, inspired by Stroock- 
Varadhan's support theorem, the following controlled ODE (from now on, we adopt the 
convention of summation over repeated indices). For each x £ T d , u £ i 1 2 oc (R+; R m ), 
let {z^(t), t > 0} denote the solution of 



dzi 1 

^) = H z (f))-- 2 

z(0) = x. 



d ^3 

dx k 



<7kj 



(zitV + aijiztfMjit), 



(2.5) 



We shall also need the controlled ODE with b replaced by b + ec, namely we shall 
denote by {z x,e (t), t > 0} the solution of 



dz- 1 
-^(t) = [bi + sci](z(fj) - - 

z(0) = x. 



dx k 



(2.6) 



We will throughout this paper make the following assumptions on the drift and the 
diffusion coefficient 



Assumption H.l The functions a, b, and c are of class C°° and periodic of period 1 in 
each direction. 

Consider now b and the o\, as vector fields on the torus T d . We say that the 
strong Hormander condition holds at some point x E T d if the Lie algebra gen- 
erated by {aj}™ =1 spans the whole tangent space of T d at x. We furthermore say 
that the parabolic Hormander condition holds at x if the Lie algebra generated by 
{dt + b} U {(TjYjLi spans the whole tangent space of R x T d at (0, x). 

Assumption H.2 There exists a non-empty, open and connected subset U of T d on 
which the strong Hormander condition holds. Furthermore, there exists to > and e 
such that, for all x £ T d , < e < £q, one has 

inf {\\u\\ L 2- zZ' E (t ) e U} < oo . 

ueL 2 (0,t o ;R m ) 

Note that, by upper semicontinuity, the supremum over x £ T d and < e < eo of the 
above infimum is bounded by a universal constant K. 

Whenever X is a random variable and A an event, we shall use the notation 



E(X; A) = E(X XA ) = E(X | A)P(A). 
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Assumption H.3 One has 

inf supE(\J?\;{r v >t})<l, 
t>o x eT d 

where V denotes the subset ofT d where the parabolic Hormander condition holds and 
Ty is the first hitting time ofV by the process {Xf}. 

Remark 2.1 One simple criteria for this assumption to hold is the existence of a time 
t 2 such that P(t v <t 2 ) = l for all x e T d . 

Denote by p e (t; x, A) the transition probabilities of the T d -valued Markov process 
{Xl' x }. We shall write p(t; x, A) for p (t; x, A). 

Lemma 2.2 Under the assumptions \H.l\ and \H.2\ the following Doeblin condition is 
satisfied : there exists t\ > 0, < E\ < £q, j3 > and v a probability measure 
on T d which is absolutely continuous with respect to Lebesgue's measure, s.t. for all 
< e < Ei, x G T d , A Borel subset ofT d , 

p £ (t 1 ;x,A)>MA). (2.7) 

Proof. It is well-known since the works of Malliavin, Bismut, Stroock et al IIMal78l 
IBis81l IKS84I IKS85I INor861 that E3] and |H72] imply that the transition probabilities 
starting from x G U have a C°° density with respect to the Lebesgue measure, and that 
this density is also smooth in the initial condition and in e. In particular this implies 
that, for every t > 0, (e, x) i— > p e (t; x, ■ ) is continuous from [0, 1] x U into the set 
of probability measures on T d , equipped with the total variation distance. Let now 
.To € U be arbitrary. It follows that there exists an open neighbourhood Uq C U of 
x, < £\ < Eq and a probability measure v such that p £ (t;x,-) > \v for every 
< £ < Ei, x G C/q- Since the measures p E (t; x, ■ ) are absolutely continuous, v must 
also be a.c. with respect to Lebesgue's measure. 

Since Uq C U and since Assumption |H.2| implies that (12. 5\ is locally controllable 
in U, it follows from the support theorem that, for every s > to and every x G T d , 
< £ < Ei, one has p £ (s; x, Uo) > 0. Since the Markov semigroup generated by the 
solutions of ( 12.3b is Feller and since Uq is open, the function (e, x) ^ p e (s; x, Uo) is 
lower semicontinuous and therefore attains its lower bound /?' on [0, ei] x T d . The 
claim follows by taking (3 = /3'/2 and t\ = s + 1, □ 

In the sequel we shall denote by {"Pet, t > 0} the Markov semigroup associated 
to the T d -valued diffusion process {Xp x , t > 0}. We shall write Vt for Vo,t- 

Corollary 2.3 Under the assumptions [77771 and \fL2\ for each < £ < £i, {X^} 
possesses a unique invariant probability measure \i e , which is absolutely continuous 
with respect to Lebesgue 's measure. Moreover, there exist constants C and q > such 
that 

sup \\p e (t;x,-)- h £ \\tv < Ce~ et , (2.8) 

x£T d , 0<£<£l 

for every t > 0. 

Proof. It follows from ( T2~7i i that for all < e < E\ 

\\V e , tl fii - P e ,tiA*a||TV < (1 - - M2||tv (2.9) 
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for any two probability measure fj,i and \x<x on T . This immediately implies the exis- 
tence and uniqueness of an invariant probability measure \i E for Xf. 

To see that fj, e is absolutely continuous with respect to Lebesgue's measure, note 
that one can decompose the transition semigroup V e ,t for t = ti as 

P e , tl =l3v + (l-0)V e (2.10) 

for some Markov operator V £ . It follows from the fact that the solutions of ( 12.21 i gener- 
ate flows of diffeomorphisms |Kun90'.|, that V™v is absolutely continuous with respect 
to the Lebesgue measure for every n. Fix now an arbitrary set A with Lebesgue mea- 
sure 0. It follows from the invariance of \i e that /i e = (V* tl ) n p e for every n > 0. 
From this, (12.101 . and the absolute continuity of v, it follows that fM e (A) < (1 — (3) n 
for every n, and therefore that fj, £ (A) = 0. 

Finally ( I2.81 l follows from iterating (|2.91 l with /.12 = Me- d 



We shall need the 



Lemma 2.4 Denote ^0 by fj» As e — ► 0, 

/i e — > /i weakly. 

Proof. The tightness is obvious, since T d is compact. Hence from any sequence e„ 
converging to 0, we can extract a further subsequence, which we still denote by {e n }, 
such that 

Me„ -> A- 

Now if / G C(T d ), t > 0, clearly V En ,tf(x) — > Vtf(x) as n — > 00, uniformly for 
cc S T d , hence 



V Sn , t f(x)fi en (dx) -» / Vtf(xMdx). 
But the left hand side equals 

f(x)fJ. Sn (dx) -> / f(x)fl(dx). 

Consequently /i is invariant under 7^* for all f > 0, hence jl = fj,, and — > /i weakly, 
as e — > 0. □ 



We can now deduce from (12.8b . using also Lemma l2~4l exactly as in Proposition 2.4 of 
|Par99l . the 

Corollary 2.5 W/ie«ever / e L°°(T d ), for any t > 0, 

in probability, as e — > 0. 
We finally assume that 

Assumption H.4 T/ie c/n/f 6 satisfies the centering condition J Td b(x)fi(dx) = 0. 
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Denoting by 

1 d 2 d 

L= 2 atj(x) dx~^ + bl(x) dx~ 

with a(x) = (T(j*(x) the infinitesimal generator of the T d -valued diffusion process 
{Xf , t > 0}, it follows from Lemma [231 below that under assumptions IH.lljH.41 
there exists a unique C 1 (T d ; R d ) solution of the Poisson equation 

Lb(x) + b(x) = , x e T d . 
This solution is given by MPV051 

poo 

b{x) = / V t b(x)dt . 
Jo 

It is not clear a priori that b is differentiable, however the following result shows 
that since b belongs to C 1 (T <i ), so does b. 

Lemma 2.6 Under assumption \H.l\ the semigroup Vt generated by ( 12.51 ) mapsC l {T d ) 
into itself. If furthermore \FL2\ and \JL3\ hold. then there exist positive constants C and 
7 such that 

"''iCHT 1 ') ^ ue '"ll/IIC 1 (T d ) 

for every f £ C 1 {T d ) such that J f{x)n{dx) = 0, and for every t > 0. 



|Pt/||cKT<*)<Ce-T*||/|| cl(T<J) , (2.11) 



Proof. The proof relies on the techniques developed in OHM0 61. Suppose that we can 
find a time <2 such that there exists a constant C and 5 > with 

\\DV t J\\ L ~ < C||/|| L « + (1 - 5)\\Df\\ L - , (2.12) 

for every test function / e C 1 (T <i ). It follows from 1HM061 Section 2] that this, to- 
gether with the Doeblin condition given by Lemma I2T21 implies that ( 12.111 ) holds. 
It remains to show that d2.12t does indeed hold. The set V can be defined as 

V = [x; Q{x) := ]T Y(x) ® Y(x) > o}, 

where y denotes the set of vector fields consisting of the aj 's, their brackets and their 
brackets with b, the number of those being bounded by d. The inequality Q(x) > 
means that the quadratic form Q(x) is strictly positive definite. We let V k be the set of 
points where Q(x) > kr l I. 

We now want to take advantage of the fact that as soon as the process X x hits 14, 
its Malliavin matrix becomes invertible. However, we face the problem that the hitting 
time of Vk doesn't necessarily have a Malliavin derivative. This is the motivation for 
the next construction. Let and f2 be two independent copies of the d-dimensional 
Wiener space, and denote by B and B the corresponding canonical processes. We first 
consider the auxiliary SDE 

dy = b(y) dt + J2 (7 j(v) dB - 1/(0) = a? , 

j 

and define as the hitting time of Vk by the process y. Let ti > be such that 

supE(|J t :c 2 |;{r^>i 2 })<l, (2.13) 

xf£T d 
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and 

-Ft*, whenf^<t 2 ; 



-oo, otherwise. 



We now define a process W by 

W t -- 



B t _ for t<T%, 

B T *+B t _ T * fbrt>7f. 

Since r^f is a stopping time, W is again a d-dimensional Wiener process. We can (and 
will between now and the end of this proof) therefore consider the process Xf to be 
driven by the process W that was just constructed. 
Consider (see ( 12.41 i) the d x d-matrix valued SDE 

TCI 

dJ? = Db(X?)J? dt + J2( D(J j( x t)Jt) dWl , Jo = I . 
We next define the (random) linear map A x : L 2 ([0, 1]; R m ) -> R d by 
A x u = J2 <•/-') A; )//,(> - 7f) , 

3=1 ■' T fc 

with the convention that = on the set {t^ = +oo}. Define the (random) map 
C x : R d — > R d by = A^A*. It then follows from |Nua95ll and the strong Markov 
property that 

Proposition 2.7 For every x G T d , the matrix C x is almost surely invertible on the set 
i T k < +°°}- Furthermore, for every k and every p there exists a constant Kk, p such 
thatE(\\C-P\\\T% < oo) < K k)P . 

Fix now an arbitrary vector £ G R d and define a stochastic process v : [0, 1] — > R' n by 

« = ^C"^ . 

(We set as before v = on the set {r^ = +oo}.) The reason for performing this 
construction is that we are now going to use Malliavin calculus on the probability space 
S = fi X equipped with the canonical Gaussian subspace generated by the Wiener 
process B. It is then straightforward to check [Nua95j that one has v G D 12 . Note 
however that since the stopping time rf is not Malliavin differentiable in general, it is 
not necessarily true that v G D 1 ' 2 if we equip S with the Gaussian structure inherited 
from W. 

This construction yields the following equality, valid for every / G C 1 (T d ) and for 
t = t 2 + 1: 

I u 11 r fc — 00 ' 

where T> stands for the Malliavin derivative with respect to B. It then follows from the 
integration by parts formula |Nua95, Sect. 1.3] that 

r 1 

m>tf{x)t = K(f(Xf) J v(s)dB( S yj - E((Df)(X?)J?£ ; r% = oo) , 
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where the integral is to be understood in the Skorokhod sense. It follows from ||Nua95l 
p. 39] that 

||^ 2 /||oo<!l^/||ooE(|Jf 2 |;K = oo}) 

+ ||/|| 00 e(/ \\v s \\ 2 ds+ [ [ \\V s v r \\ 2 drds 
vo Jo Jo 

Since, for every k, the second term is bounded uniformly in x 6 T d , it remains to show 
that we can choose k such that 

supE(|Jf 2 |;{7f = oo})<1. 

X 

Assume this is not the case. Then to each k, we can associate a point x& such that 

E(|Jf a *l;K*>t a })>i. 

We can assume (up to extracting a subsequence) that Xk — > x. For all £ G N, 

1 < lim supE(|J**|; {t** > t 2 }) < lim supE(| J%\; {rf > t a }) 

<E(|JfJ;K>t 2 }), 

since the mapping x — > X{r| : >t,} is a - s - upper semicontinuous, as the indicator func- 
tion of a closed subset of trajectories. Since 

/_™X{Tf>t 2 } = X{r»>i 2 } a-s., 

this contradicts (12.131 1. □ 

Remark 2.8 We claim that the assumption lH.3l is close to being sharp for the conclu- 
sion of Lemma l2~6l to hold. The simplest example (which is however not a diffusion 
on the torus) on which this can be seen is as follows. Let r be an exponential random 
variable with parameter b > 0, and 



X? 



e at x, if t < r; 
0, otherwise. 



In this case 

E(|Jf|; {r > i})= e (a - &)t , 

and 

||^/!loo=e (a - fc)t ||£/||oo, 
so that ||Pf/|| c i — > if and only ifJOMds. 

3 The homogenization result 

The goal of this section is to show the 

Theorem 3.1 We have that, in the sense of weak convergence on the space C(R+) 
equipped with the topology of uniform convergence on compact sets, 

r^I 1 , where Xf = x + Ct + A 1/2 W t , 
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as e — > 0. //ere, {W^t, t > 0} is a standard Brownian motion, and the homogenized 
coefficients C and A are given by 



C= / (I + V6)c(x) n{dx) , 

A= [ (I + Vb)a(I + Vb)*(x) n(dx) 



In order to prove this theorem, we need to get rid of the term of order e 1 in the SDE 
(12. 11 1. The trick is to replace X^' e by 



X? > £ := Xf< e + eb 



X^ e 



and to take advantage of the 

Lemma 3.2 The following equality holds almost surely: 



X x t ' e =x + eb(~) + J V +V ® C (^J~) ds + J (I + ^ b)<7 {^~) dWs 
Proof. Let q : R d — ► R + be a smooth function with compact support, such that 

q(x) dx — 1, 

and define g rl (x) = n d g(nx). We regularize 6 by convolution 

b n = b* Q n - 

We now deduce from Ito's formula that 



<re,n v x,e , i / 

x t -= x t +s h 



e 



x + sbn^-j +e- L J (Lb n + b)\^—jds + J (I + Vb n )c\^-jds 

(I + V6„W -J-J dW a . 

We want next to let n — > oo. Clearly 6„ — > 6 and V6„ — > V6 pointwise, and the two 
sequences are bounded, uniformly with respect to n and x E T d . It remains to treat the 
term containing the second order derivatives. One has 

Lb n = (Lb) * g n + <p n , 

and since b is a weak solution of the Poisson equation, 

(Lb) * Q„ = -b * Q n -» -b, 

the sequence being again uniformly bounded. It remains to study the sequence ip n . 
Using again the convention of summation over repeated indices, we have 



If d 2 b 



f n (x) = 2j [oij(%) - aij(x - y^ Q x .Q x . ( x ~ v^ n Q( n y) d y 

+ [ lh(x) - h(x - y)]^-(x - y)n d g(ny)dy. 

JR d OXi 



(3.1) 
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The fact that the second integral in ( 13. U converges to and is uniformly bounded is 
easily established. We now integrate by parts the first integral, yielding 

d 2 b 

[a>ij(x) - a i3 (x - y)]- — —(a; - y)n d g(ny)dy 
R d oxiOXj 

dctij db j 

-7; — (x~y)- — (x-y)n g(ny)dy 

R d OXi OXj 

+ [ [ai 3 (x) - aijix - y)]-^-(x - y)n d+1 g' l (ny)dy 
J R d dXj 

The first term on the right hand side is uniformly bounded and converges to 

daij db 
-(x)—{x) . 



dxi dxj 



The second term is equal to 



/ V ■ Va i:) (x - y')^-{x - y)n d+1 g' l (ny)dy , 
J Rd dxj 

where \y'\ < \y\. This last quantity is equal to a bounded sequence converging to 0, 
plus 

da^ , „ db , , f rf , 1 , x , dan , ■. db , ^ 
-^l(x)—(x) / y k n d+1 g' l {ny)dy = -—^{x)—{x). 



dxk dxj J dxi dxj 

The lemma is established. □ 

We now proceed with the 

Proof of Theorem \3.1\ We first note that since 

< Ce, 

the Theorem will follow if we prove that, as e — > 0, 

X x ' e ^X x . 

But since 



X^ E = x + J\l + Vh)c(^\ d 



s + Ml 



one has 

((M e ))t = J (I + V6)o(J + V6)* f ^) ds . 

and the result follows from Corollary 12. 5| and the martingale central limit theorem, see 
e.g. EK86l Thm 7.1.41. □ 



We conclude this section with the following result, which extends the averaging result 
of Corollary [33] It is going to be needed in the applications of sections 5 and 6. 
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Proposition 3.3 Let f G C(R d x T d ). Then for every t > 0, the following convergence 
holds in law as e — ► : 



where 



fix) 



ds 



f 

Jo 



f(x,y) Kdy)- 



Proof. It is easily checked that / is continuous, hence 



f(X x s )ds, 



as e — > 0. It then suffices to show that 



s •> 



X?> E 



f(X? E )ds 







in probability, as e — > 0. This is proved by exploiting the tightness in C(R+) of the 
collection of processes {X x ' £ , e > 0}, and Corollary 12. 5 1 For the details, see the proof 
of Lemma 4.2 in MPar99ll (there is a misprint in the statement of that Lemma). □ 



4 Characterization of the image of the homogenized diffusion ma- 
trix 

The aim of this section is to give a characterisation of the range of A which can very 
easily be read off the support of the invariant measure fj,. We first state all of our results 
and then prove them in the following subsection. 

4.1 Statements of the results 

It will be convenient in our statements to choose an arbitrary particular point xq g U. 
We define the set T, consisting of loops on T d , starting and ending at xq, which take 
the form 

7(a) = <°(s), < s < t, 

with arbitrary t > and u G L 2 (0, t; R m ). We call loops in T admissible. It follows 
from IH.21 that the set of admissible loops does not depend on the particular choice of 

X . 

We next define a mapping g : T — > Z d as follows. Lifting each loop 7, so as to 
transform it into a curve 7 in R d , we define 

5(7) = lH) - 7(0), 

where t is the end time of the loop 7. The first step in our characterization of the image 
of A is to show that 

Theorem 4.1 Under the assumptions \H.l^L4\ g(T) C ImA. 
The converse to this result takes the following form: 
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Theorem 4.2 Under the assumptions \[Ll\\H.4\ for each e € ImA, with \e\ = 1, and 
for every S > 0, there exists a loop 7 £ T such that 



3(7) 



< 5. 



|S(7)| 

It follows immediately from the above theorems that under the same assumptions, 
Corollary 4.3 ImA = span{g(7), 7 e T}. 

Corollary 4.4 The matrix A is non degenerate if and only if there exists a collection 
of admissible loops {71, . . . , 7d} C T such that 

span{ 5 (7i),...,<7(7 rf )} =R d . 

The next corollary is slightly less obvious and will be shown in the next subsection. It 
will mainly be useful in the proof of Theorem l4.6l below. but since it is not a priori an 
obvious fact, we state it separately. 

Corollary 4.5 The set G = {3(7), 7 € T} is a subgroup ofZ d . 

We denote by S the interior of the support of the invariant measure fi. It follows 
from IH.21 that the support of /i is then the closure of S. It turns out that as far as the 
characterisation of ImA is concerned, we can replace admissible loops by arbitrary 
loops in S. This may sound somewhat surprising at first sight, since it is certainly not 
true in general that every loop in S is admissible. One has however 

Theorem 4.6 Every admissible loop 7 satisfies 7 C S. Conversely, for every loop 
Q C S, one can find an admissible loop 7 with g(g) = 5(7). 



Remark 4.7 Combining Theorem 14.61 with Corollary 14.31 we find that 1mA is com- 
pletely characterised by the topology of the set S. This is however not necessarily 
true if we drop the regularity conditions IH.2I - IH31 In particular, it follows from Corol- 
lary [43] that one can find a basis for ImA which has rational coordinates. A 'coun- 
terexample' (with d = 2 and m = 1) to this claim is given by taking b = c = and 
a 1 = (1, v2), so that ImA is the vector space generated by (1, v2). Of course, this 
example satisfies neither lH31 nor llL3l and does therefore not contradict our results. 

Remark 4.8 The argument on page 264 in [IK094| shows that A is non degenerate 
whenever the set span{ .9(7), 7 £ So} is equal to R d , where Sq is a connected subset of 
the support of the invariant measure \i where a is elliptic. Clearly the result presented 
here is stronger. 

4.2 Proofs of the results 

We start with the the 

Proof of TheoremWJ] Let £ € R d be such that (A£,£) = 0. We will prove that for 
any 7 S T, (3(7), £) = 0. We shall make use of the notation b^(x) = (b(x), £), and of 
the three following facts : 

• 65 solves the Poisson equation 

Lb 6 (x) + & e (a;) = 0, xeT d ; (4.1) 
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(A£, £) =0, which, when we explicit the matrix A, amounts to the fact that 

(Vb € ,a.j) + (a.j, £} = 0, 1 < j < d, for ^ a. e. x. (4.2) 



Note that ( 14.21 ) holds true /i-almost everywhere. It can be noted that the trajec- 
tories of the solutions of ( 12.51 ) starting from the point xq remain in the interior 
of the support of /i, which is absolutely continuous with respect to the Lebesgue 
measure. Hence ( 14.2b holds true almost everywhere in the neighbourhood of 
such a trajectory, consequently everywhere there, due to continuity. 

• We can now differentiate ( 14.2b with respect to Xk, and multiply the resulting 
identity by a k j , from which we deduce that 



a k j(x) 



d 
dxh 



db t 

' dxi 



0*0 = a k j 



doj. 
dxi 



(x), x G T d . 



(4.3) 



Choose a loop 7 G T, and denote by t its end time. We now have, for < s < t, 



= <6,V6 c )(7(«))-5 



-6«(7(s)) 



as 



db^ d<Jij 
dxi dxk 



(700) + 

<9(7y Sfcj 



06* 



(7(s))uj(s) 



13 dxidxj ' 9a;fe 9xj 



(7(s)) - ^ij(j(s))Uj(s) 



d(Ti 



(Jkj 



dx k 



(7(«)) - £i0y(7(s))Mj(s) 



where we have used ( 14.1b and (14.2b for the third equality, and ( 14.3b for the fourth equal- 
ity. Integrating from s = to s = t, we deduce that 



(ff(7),0 + (^(7(i))-^(7(0)),0 =0, 
from which it follows, since 6 is periodic, that (5(7), £) = 0. The result follows. 



□ 



Proof of Theorem \4.2\ Let e 6 ImA with |e| = 1 and a > 0. We denote by B(x,a) 
the ball in R d of radius a, centered at x. Note that 

P(^ 1/2 Wi G B(e,a)) > 0. 

It follows that for e > small enough (here X 6 is defined with c = 0), 

Consequently 

It then follows from Stroock-Varadhan's support theorem that there exists a control 
u G i 2 (0, 1/e 2 ; R m ) such that the corresponding z-trajectory, lifted to R d , satisfies 



1(0) = i , \z(l/e 2 )-(x + e/e)\<a/e. 
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Now from lH.2l there exist t < t and u £ L 2 (-±r, ^ + t;R m ), such that z{j? + t) = 
x o + <?(7)> where 

7 = |z(s), 0<s<-^+i}er. 
We have constructed a loop 7 such that 



e 

3(7) - - 

e 



e 



for some universal constant C (see the remark following assumption lH.2l ). Multiplying 

by e, we get 

\eg(l) ~ e| < a + eC, 

and consequently 

g(7) 
15(7)1 ' 

from which the result follows, provided we let e < 8 /4C and choose a = 8/4. □ 



< 2(a + eC), 



Proof of Corollarv \4.5\ It follows from Theorem 14.21 that linear combinations of ele- 
ments of G with positive coefficients are dense in ImA Since they form a cone, they 
must therefore coincide with ImA Denoting the dimension of ImA by k, it follows 
from Theorem 14. II that, for any g £ G, there exist ai, . . . ,a.k £ R+ and linearly 
independent elements gi,. ■■ ,gk £ G such that 

k 

-.9 = 

i=i 

In fact, since g and all gi have integer coordinates and since the inverse of a matrix with 
integer coefficients has rational coefficients, the a/s must be rational. Consequently, 
there exists n £ N such that — ng £ G, from which it follows immediately that —g £ 
G. We have proved that G is a group. □ 



Proof of Theorem \4.6\ The first statement follows from the fact that any admissible 
trajectory starting at xq stays in S. 

Conversely, we first note that since S is connected, to each loop g C S, we can 
associate a loop g' C S starting from xq, with g(g') = gig). So we might as well 
assume that g starts from xq. Denote by S the lift of S to R d . It now follows that xo 
and xq + g{g) belong to the same connected component S' of S. It remains to show that 
there exists an admissible loop 7 such that 7 joins xq and xq + gig). In other words, it 
remains to show that gig) £ G. 

For that sake, first notice that 

S = U/jgzdCCfc), 

where C(fc) denotes the set of accessible points from xq + k by the solution of the 
controlled ODE ( 12.51 ) lifted to R d . Consequently we have that 

S = ^>kez d , xo+keS'Cik). 

Recall that xq + gig) £ S 1 . Since each C(fc) is open, there exists n > and a chain 
/1 ,...,/„ of points in Z d such that, with fo = and f n +\ = gio), 



Cifi) n Cif l+1 ) + < i < n. 
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Pick a point y £ C(/j) n C(/j+i). There exist a control such that the solution of 
(12.5b (lifted to R d ) starting from xq + /; reaches y in finite time, and another one 
such that the solution starting from xq + also reaches y in finite time. It follows 
from assumption IH.2I that there exists a control such that the solution starting from y 
reaches some xq + f in finite time. Hence / — fa, f — /j+i £ G. Since G is a group, 
fi+i - fi G G, for all < i < n. We have proved that g(g) G G. □ 

5 Homogenization of an elliptic PDE 

In this section, we show how to apply the homogenization results obtained in the previ- 
ous sections to the elliptic homogenization problem ( II. lb . Let D be a bounded domain 
in R d with a C 1 boundary, and define 



We assume that / £ C(R ), that it is periodic with respect to its second variable, and 
that there exists 8 > such that 



Remark 5.1 We will need conditions ( 15. Il l and ( 15.21 i in order to deduce some uniform 
integrability. Condition ( 15. U must be checked in each particular example. It is always 
satisfied with a = (unless one is in the case A = 0), in which case (15.2b requires that 
/(x, y) < 0, x e R d , y g T d . This condition can be relaxed only if (15.11 ) is satisfied 
with some a > 0. 

The solution of the elliptic PDE 



t £ = M{t > 0, Xf D}. 



Let a > be such that for all x £ D, 



supEa; exp(ar e ) < oo. 

£>0 



(5.1) 



/(a:,y)<(a-<5)+, Wx £ R d , v € T d . 



(5.2) 




is then given by the Feynman-Kac formula 




We define as before 






dxidxj 
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From Theorem l3.ll X s => X as e — > 0, where 

X t = x + Ct + A 1/2 W tl t > 0. 
We assume here moreover that 

A is strictly positive definite. (5.3) 

Then 

t = mf{t > 0, X t <£ D] 

is an a.s. continuous function of the limiting trajectory {Xt}. It follows from Proposi- 
tion |33] that 



f\XI,-±) ds^ J D(X s )ds. 

From (15. It et (15. 2t we deduce the necessary uniform integrability in order to establish 
the 

Theorem 5.2 Under the conditions\Hl\ BO] gO] El <ED, <E3 ««<^ <E3, 



u e (x) 



D(X S ) ds 



, as e — > 0, 



where u{x) := E x \g{X T ) exp (J^ D(X S ) (is)] i's f/ie solution of the elliptic PDE 

j Cu(x) + D(x)u(x) = 0, x e L>; 

u(s) = <?(a;), a; £ <9Z), 

af Zeasf ;« f/ie viscosity sense. 

Remark 5.3 The assumption that A be nondegenerate is not really necessary for our 
purpose. All we need is that for almost all (with respect to the law of X!j!) points 
x G dD, (An(x),n(x)) > 0, where n(x) denotes the normal at x to dD. Depending 
on the geometry of D, this can be true even if dim(Imj4) = 1 and d > 1. 



6 Homogenization of a parabolic PDE 

Assume that a, b, c satisfy the above assumptions. Let e € C 1 (R d , R), / € Cb(R d x 
R d , R), e is periodic, / is periodic with respect to its second argument, and g € C(R d ) 
grows at most polynomially at infinity. For each e > 0, we consider the PDE: 



(t, x) = C E u £ (t, x)+(-e(-) + f(x, -)) u £ (t, x) 



We assume that 



Define 



u £ (0,x) = g(x),xe R d . 

e(x)p(dx) = 0. 



(6.1) 



(6.2) 



Yt 



i fx 



ds. 
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Then the solution of d6.U is given by 

u"(t,x)=E[g(X[)cxp(Y t £ )] , 

where Xf is the solution of the SDE (12. 11 1. 
Define 

poo 

e(x)= / E x [e(X t )]dt, 
Jo 

poo 

k(x)= / E x [bi(X t )]dt, 1 = 1,.. 
Jo 

the weak sense solutions of the Poisson equations 

Le(x) + e(x) = 0, 
Lb l {x) + bi(x) = 0, i=l,...,d. 



(6.3) 



We let 



A= (1 + Vb)a(I + Vb)*(x)n(dx) 

Jjd 

C = I (I + V&)(c + aX7e)(x)fi(dx) 

D(x) = I (ive*aVe+/(a;, ■) + \7ec)(y)n(dy) 
Jjd z 

Then, defining again X? = x + Ct + A*W t , 

u(t,x) = E[g(X?)eti D{X *' )ds ] 
is the solution (at least in the viscosity sense) of the parabolic PDE 
1 



(6.4) 



dt (t,X) ~ 2 Al3 d Xi dx 
w(0, x) = g(x), x G R 



(t, x) + C l — {t, x) + D(x)u(t, x), 
ox. 



Theorem 6.1 Under the conditions \H.1\ \H.2\ \H.3\ \H.4\ the smoothness of e, the 
boundedness and continuity of f, the continuity and growth condition on g and ( 16.21 ), 
for aIlt>0,x<E R d , 

u e (t, x) — > u(t, x) 

as e — » 0. 

For the proof of Theorem l6.ll we need a result which is proved exactly as Lemma l3~2l 
namely 



Lemma 6.2 Define 



y/ = y/ + £ e(^). 



Then the following holds for each e > and t > : 



*f -(f) 



[/ + Vec] X 



X 



ds 



Via 



dW* 
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We can now proceed with the 

Proof of Theorem \6.1\ We first define a new probability P by the formula 



dP 
dP 



= exp 



A'f 



1 



Veer — s - )dW s - - / (Ve, aVe) — ) ds 



We next remark that it follows from (16.3fr that the behaviour as e — > of u s (t, x) is the 
same as that of 

u 6 (t,x) := E[ 5 (l t e )exp(Y/ - ee(x/e))]. 
The definition of P then yields 



u e (t, x)=E 



g(X!)exp I f[X e a) ^)ds 



where we defined 



/ = / + Vec+-(Ve,aVe). 
On the other hand, it follows from Girsanov's theorem that 



Xf=x- 



e&Q + J (I + V&)(c + aVe) f^-\ ds 



( £L ) 



where Wt = Wt — J * cr* Ve ^^^) ds is a Brownian motion under ] 



An obvious adaptation of the proof of Theorem l3 . 1 1 shows that X x ' e X x , where 
Xf — x + Ct + A 1 / 2 W t , C is given by i6.4\ , and {Wt} is a Brownian motion under 
P. Since it follows from Proposition |3.3| that 



f[xt 



X! 



ds 



WXf) ds 



as e — > 0, the result follows. 



□ 



7 Examples and counterexamples 

This section provides several examples of degenerate diffusions for which our results 
apply. In each case, we furthermore give a characterisation of the range of the effective 
diffusion matrix A. 

The figures should be interpreted as follows. The little black arrows show the vector 
field b. The shaded dark grey regions denote the points where the strong Hormander 
condition holds. 

Our first example (which is the only one that we are going to work out in some 
detail) is a typical example of the type of diffusions for which our conditions IH. H - IIOl 
apply. Observe first that the diffusion {X t , t > 0} on T d has the measure with density 
p as invariant measure, and lH.4l is satisfied if and only if there exists a mapping H from 
T d into the set of d x d antisymmetric matrices, such that 

bi(x) = ^— ir-(paij + Hij){x\ 1 < i < d, x e T d . 
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Let a : T d — » [0, 1] be a smooth function such that the set A = {x \ a(x) = 0} 
has a finite number of bounded connected components. For g > 0, we define the set 
A g := {x, d(x, A) < g} and we assume that there exists g > such that A g intersects 
neither {xi = 0} nor {xi = 1/2}. 

We choose our diffusion matrix a to be given by a(x) = a(x)I, and we let 



/ 



H = 





\-(27r)" 1 cos(27ra;i) 



Hence, if we choose 



2b(x) = V(a 2 )(x) + 



(27r)- 1 cos(27T2;i)\ 




( ° \ 



\six\(2iTXi)) 



then the Lebesgue measure is invariant for the process Xf . The following picture 
shows the vector field given by b, together with an example of a set A (the white disks) 
that satisfies our conditions. 



I 



li 



Assumptions IH.ll and IH.4I are satisfied by construction. Assumption IH.2I is satisfied 
since, by taking the trivial control u = 0, there exists a time t such that every solution to 
the controlled system reaches the complement of A before time t. The same argument 
shows that lH.3l is satisfied as well. 

Our second example is one where the diffusion matrix of the limiting Brownian 
motion is non-degenerate even though the area in which the strong Hormander condi- 
tion holds does not intersect the boundaries of the fundamental domain. In particular, 



Examples and counterexamples 



20 



the argument of |JKQ94| mentioned in Remark |4~81 does not cover this situation: 




In this case, it is easy to see that every point of the torus can be reached by the diffusion 
on the torus, so that \i has full support and therefore A has full rank. 

Our third example is in a way the opposite of the first one. It shows that it is possible 
for the limiting Brownian motion to have zero diffusion coefficient, even though the 
area in which the strong Hormander condition holds stretches over the whole space 



/ ' » \ 

/ - - \ 


/ ' *■ \ 
/ ■ - \ 


t ' > \ 


li 


\ - - / 

s v / 

/ x x \ 


\ - - / 
\ \* / 

/ X X \ 


\ - - /• 

\ X X / 




t<*\ 

/ • - \ 


/ v\ \ 
✓ \ 


/ X X ^ 

x* • - \ 


✓ -\ 


\ - - / 

1 v/ 

/ X X \ 


\ \' ; 

/ x X \ 


\ 

\ w / 

/ X X \ 


\ — ✓ , 
\ w / 



In this case, the support S of the invariant measure is given by the gray disk, so that 
5(7) = for every loop in S. 

We finally give an example where the effective diffusivity degenerates in one direc- 
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tion. 




In this particular case, the range of A is the one-dimensional subspace of R 2 spanned 
by the vector (1, 2). 
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